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ABSTRACT 

We discuss the observations and theory of star cluster formation to argue that clusters form dy- 
namically cool (subvirial) and with substructure. We then perform an ensemble of simulations of 
cool, clumpy (fractal) clusters and show that they often dynamically mass segregate on timescales far 
shorter than expected from simple models. The mass segregation comes about through the produc- 
tion of a short-lived, but very dense core. This shows that in clusters like the Orion Nebula Cluster 
the stars > 4Mocan dynamically mass segregate within the current age of the cluster. Therefore, 
the observed mass segregation in apparently dynamically young clusters need not be primordial, but 
could be the result of rapid and violent early dynamical evolution. 

Subject headings: methods: N-body simulations — stars: formation — stellar dynamics — galaxies: 
star clusters 
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1. INTRODUCTION 

Most stars are observed to form in clusters (proba- 
bly 70—90%; Lada & Lada 2003.) Therefore, under- 
standing cluster formation is the key to understanding 
most star formation. Unfortunately, observations of very 
young clusters (<IC 1 Myr) are hampered by extinction 
due to those clusters still being embedded in their na- 
tal molecular clouds. Most often we are restricted to 
observing older clusters which is problematic as these 
clusters may well have undergone significant dynamical 
evolution. In particular, young clu sters seem to expan d 
significantly in the first few Myr (jBastian et al.l f2008f ) . 
possibly due to the effects of multiple stellar encounters 
()van den Berk ct al. 20 03) or by the rapi d expulsion of 
primordial gas (Good win fc Bastianll2006f l. 

A potentially significant observation is that many 
young clusters appear 'mass segregated,' i.e., the most 
massive stars are concentrated towards the center of the 
cluster (Hillenbrand & Hartmann 1998; de Grijs et al. 
2002a, b,c; Gouliermis et al. 2004). Clearly, mass segre- 
gation could be primordial - clusters may form with the 
most massive stars concentrated at or near the center. 
Alternatively, mass segregation could be dynamical - the 
most massive stars migrate into the center of the cluster 
after formation due to two-body interactions. As well as 
potentially providing constraints on clustered star forma- 
tion, the origin of mass segregation may help distinguish 
models of massive star formation. In particular, are the 
masses of the most massive st ars set by the mass of the 
core in which they form (e.g.. iKrumholz et al.ll200"7l ). or 
by competitively accreting mass due to a favorable po- 
sition in the cluster (e.g., Bonnell et al. 1998; see also 
Krumholz et al. 2005, and Bonnell & Bate 2006)? If 
mass segregation is primordial it might well argue in fa- 
vor of competitive accretion, as in that scenario massive 
st ars should form in the ce nter of a cluster. 

iBonnell fc Davie^ ()1998l ) investigated the evolution of 
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virialized and subvirial (thus collapsing) smooth Plum- 
mer spheres. They showed that dynamical mass segre- 
gation cannot occur rapidly enough in these situations 
to explain the observations of mass segregation in young 
(few Myr old) clusters such as the Orion Nebula Clus- 
ter (ONC). In this letter we will examine the dynami- 
cal evolution of initially cool and clumpy (fractal) stellar 
distributions. We will show that a combination of sub- 
virial velocities and substructure can often lead to dy- 
namical mass segregation on a timescale that is much 
shorter than would usually be expected. In § 2 we argue 
that cool, clumpy initial conditions are realistic - indeed, 
they are both observed and expected from theory. In § 3 
we demonstrate that these initial conditions can lead to 
rapid dynamical mass segregation. In § 4 we discuss the 
implications of this result and draw our conclusions. 

2. THE INITIAL CONDITIONS OF STAR CLUSTERS 

In this section we review the observations and theory of 
star cluster formation in molecular clouds. Our aim is to 
show that both observations and theory lead us to expect 
that clusters form both subvirial and with substructure. 
Star clusters are initially sub structured. Stars 
form from dense cores in molecular clouds (e.g., 
I Ward-Thompson et al.l [20071 ). Molecular clouds are ob- 
served to have significant levels of substructure in both 
density and kinematics (Williams 1999; Williams et al. 
2000; Carpenter fc Hodapp 2008). This is not surprising, 
as we believe that molecular cloud structure is dominated 
by supersonic turbulence which will produce complex 
structures (see, Mac Low fc Klessen 2004; Ballesteros- 
Paredes et al. 2007, and references therein). It is not 
surprising, therefore, that very young (< 1 Myr) star 
clusters also show significant levels of substructure (Lar- 
son 1995; Testi et al. 2000; Elmegreen 2000; Lada & 
Lada 2003; Gutermuth et al. 2005; Allen et al. 2007) 
Detailed statistical studies also show the presence 
of substructure in young clust ers, and that substruc- 
ture is ra pidl y eras ed (Cartwri ght fc WhitworthI |2004 
ISchmeia et al.il2008f ). 

Star clusters form subvirial (cool). There is increasing 
evidence that star clusters are born subvirial. Prestellar 
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cores are often observed to move subsonically (Belloche 
et al. 2001; Andre 2002; Walsh et al. 2004; Peretto et 
al. 2006; Kirk et al. 2007). In hydrodynamic simulations 
of cluster formation starting from dynamically hot (glob- 
ally unbound) gas the stars that form also appear to have 
subvirial motions (see e.g. Bate et al. 2003; Bonnell et 
al. 2003). Proszkowet al. (2009) also find that the Orion 
star forming region shows signatures of subvirial dynam- 
ics on scales of about 10 pc (see also Adams et al. 2006). 
Such subvirial motions might well be expected from a 
turbulent model of star formation in clouds where stars 
form in converging/colliding flows (see references above, 
also Adams et al. 2006). It is also worth noting that sub- 
virial initial conditions are required to erase substructure 
as rapidly as is observed (Goodwin & Whitworth 2004). 

3. THE EVOLUTION OF COOL, FRACTAL CLUSTERS 

In this section we simulate clumpy (fractal), cool clus- 
ters and show that they mass segregate dynamically in a 
time similar to the crossing time of the initial (or final) 
system. 

3.L Initial Conditions 

We simulate the dynamical evolution of star clusters 
containing 1000 single stars with masses sampled from a 
three-part power-law mass function (Kroupa 2002), with 
minimum and maximum masses of O.OSMQand SOMq, re- 
spectively. The stars initially have a fractal distribution 
within a sphere of radius 1 pc, with initial velocities such 
that nearby s tars have simila r velocities, as described in 
detail by Go odwin fc Whitwo rth (2004). We neglect the 
effects of stellar evolution because of the short duration 
of the simulations (4 Myr) . 

In this letter, we restrict our investigations to clusters 
with a fractal dimension of 1.6 (giving a very clumpy dis- 
tribution, where 3.0 gives a uniform sphere) and a virial 
ratio of Q — 0.3, where Q is the ratio of the kinetic to 
the (modulus of the) potential energy (so that Q — 1/2 
is virialized). Such initial conditions give the most ex- 
treme dynamical evolution and the most rapid dynamical 
mass segregation. We will consider a much fuller range 
of parameter space in a follow-up paper (R. J. Allison et 
al., in prep.). Here, we just wish to illustrate that rapid 
dynamical mass segregation can occur with plausible ini- 
tial conditions, and how this happens. The simulations 
were carried out using the k ira integrator in STARLAB. 
(jPortegies Zwart et al.|[2001h . 

3.2. Results 



Figures 1(a) and 1(b) show the stellar distributions 
initially, and after 1 JVIyr of dynamical evolution, re- 
spectively. A comparison of the plots clearly shows that 
the cluster has evolved from a clumpy and non-mass seg- 
regated state to one which has erased substructure and 
appears to be mass segregated. 

We apply the method of lAUison et al.l (|2009D . which 
compares the minimum spanning trees (MSTs) of high- 
mass stars to those of a random selection of stars to pro- 
duce a quantitative measure of mass segregation. If the 
MST of the N most massive stars is significantly shorter 
than that of a number of sets of N random stars then the 
cluster is mass segregated. The degree of mass segrega- 
tion can be quantified by the ratio of the lengths of the 
average randomly selected star MST to the most massive 
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Fig. 1. — Typical stellar distributions at and ~1 Myr. The 
triangles indicate the positions of the 10 most massive stars. In this 
run the most massive and tenth most massive stars have masses of 
23MQand 3.5Mq, respectively. 



Star MST, A fsee lAllison et al.ll2009l for details). The 
greater A is relative to unity, the more mass segregated 
a cluster is. 

Figure [2] shows the evolution of A for four subsets of 
the N = 10, 20, 50 and 100 most massive stars in the 
cluster. The cluster is not mass segregated initially (A = 
1), but after 1 Myr the 10 most massive stars develop a 
significant level of mass segregation (A ^ 3). The error 
bars in Figure [2] represent the instantaneous standard 
deviation at each simulation snapshot.^ In Figure [2] we 
can also see that the 20 and 50 most massive stars also 
mass segregate, but by a much smaller amount. Beyond 
the 50 most massive stars little mass segregation is seen. 

The cluster illustrated in Figure [1] is a fairly typical ex- 
ample of the evolution of cool, highly fractal clusters: a 
collapse from the cool initial state, which erases substruc- 
ture and also imprints mass segregation. Unfortunately, 

The 10 most massive stars are mass segregated at a similar level 
for more than a crossing time after the initial violent relaxation 
phase. In a crossing time the cluster can completely mix and every 
star in the cluster can migrate to any other position in the cluster. 
This means that any feature that remains constant over a crossing 
time is a real feature and if the evolution of the cluster were also 
accounted for the significance would be much greater than shown. 
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Fig. 2. — The evolution of the mass segregation ratio A for the 
10, 20, 50 and 100 most massive stars in the cluster (from top to 
bottom). The envelopes indicate \a deviations, no error is shown 
for TV = 20 and 50 for clarity. The dotted line indicates a A of 
unity, i.e., no mass segregation. 



and unavoidably, using fractal clusters introduces a cer- 
tain degree of randomness as each fractal - while formally 
the same (i.e., the same fractal dimension and virial ra- 
tio) - is very different in its initial distribution. Thus, 
when dealing with fractals, it is vital to perform a large 
ensemble of simulations. 

We simulate 50 different clusters, varying only the ran- 
dom number seed used to create the fractal and initial 
mass function. Of these 50 simulations, 29 mass segre- 
gate within 1 Myr, and 44 show mass segregation within 
4 Myr. Only 6 show no significant mass segregation by 
the end of the simulations at 4 Myr. Here, mass seg- 
regation is defined as any mass segregation event that 
lasts for longer than 0.1 Myrs and for which A has a 
significance > 1. 

3.3. Mass segregation mechanism 

The mass segregation observed in these simulations is 
from a wholly dynamical origin, and arises from the col- 
lapse of the cluster. The initial conditions used in these 
simulations place the cluster far from equilibrium. The 
cluster undergoes a gravitational collapse and violent re- 
laxation phase very early on in its evolution in an at- 
tempt to virialize itself. 

This collapse creates a dense core containing roughly 
half the mass in a radius of only ~ 0.1 pc. The dense 
core only lasts for 0.1 - 0.2 Myr, but due to its small size, 
this is around 10 - 20 crossing times of the core, making 
it dynamically old. 

Whilst 10 - 20 crossing times is not enough time to 
reach full equipartition, it is suffi cient time to m ass seg- 
regate the most massive stars. ISpitzed ()1969l ) showed 
that the mass segregation timescale for a star of mass 
M, iscg(M), is 



where (m) is the average mass of a star in the cluster 
{'^ 0.4Mo for a typical initial mass function), and treiax 
is the two-body relaxation timescale of the cluster. The 
two-body relaxation timescale depends on the number of 
stars in the cluster, N, and the crossing time, Across, as 
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where the crossing time is simply the radius of the cluster 
R divided by the average velocity of a star, a. 
Thus, Eq. [T]can be rewritten as 
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Typical values for these parameters for the dense cores 
are TV - 300 - 500, i? ~ 0.1 - 0.2 pc, (m) = O.4M0, 
and cr ~ 2 km s""'^. The lifetimes of the dense cores 
during which they can mass segregate are tscg ~ 0.1 Myr. 
This suggests that clusters will be able to mass segregate 
above M ~ 2 - 4Mq. 

In Figure [2] we show that below about the 50th most 
massive star there is no further mass segregation. In this 
simulation, the mass of the 50th most massive star is ~ 
2Mq - in good agreement with the calculation presented 
above. 

Interestingly, lAllison et all (|2009D find that the ONC 
appears to be mass segregated do wn to ^ 5Mp) but 
not below that mass. Hillenbrand fc HartmannI ()1998l ) 
also find evidence for mas s segregation around 5M0 , and 
iMoeckel fc BonnelH (|2009D find that they cannot explain 
the observations of the ONC unless only the most mas- 
sive stars are mass segregated. This is exactly the sit- 
uation predicted in our simulations - mass segregation 
down to a few solar masses, but not below. Indeed, ob- 
servations of the mass down to which a cluster is mass 
segregated (and the dynamical age of the cluster) should 
provide constraints on the density and duration of the 
dense phase undergone by the cluster and hence the ini- 
tial conditions of the cluster. 

The mass segregation seen in these simulations is due 
to the collapse of the cluster and formation of the dense 
core. Some early, low level, mass segregation is ob- 
served in the simulations, due to the dynamical evolution 
and merger of the sub -clumps in the initial distribution 
(jMcMillan et al.ll2007h . but the presence of long lived, 
high level mass segregation is due almost entirely to the 
later evolution of the dense core, not to prior mass seg- 
regation in clumps. 

The reason that clumpy subvirial clusters are able to 
mass segregate whilst smooth subvirial clusters do not 
(such as those simulated by Bonnell & Davies 1998) is 
due to clumpy clusters being able to collapse to a far 
denser state than smooth clusters. The potential energy, 
SI, of a cluster of mass Afdus and radius R is given by 



R 



where a is a structure parameter (for example in Plum- 
mer sphere if R is the Plummer radius then a ^ 0.75). 
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For a cluster with an initial potential energy Qq (with 
radius Ro and structure parameter ao), and a final poten- 
tial energy fij (with radius Rf and structure parameter 
af) the initial and final potential energies are related to 
the initial virial ratio Qq by 

f^o(i-go) = ^, 

assuming the cluster ends virialized. The ratio of initial- 
to-final radii is 

Rf aj 

In the case where a cluster starts smooth (e.g. a Plummer 
sphere) it will retain this structure and so ao ^ af and 
so if a cluster starts with a virial ratio of Qo = 0.3 it will 
only collapse by a factor of ^ 1.4. This is not enough to 
significantly increase the speed of mass segregation. This 
is consistent with the lack of mass segregation found by 
Bonnell & Davies (1998) in the collapse of a Qo = 0.1 
Plummer sphere (which would collapse by a factor of 
~1.8). 

However, when a clumpy cluster collapses it erases sub- 
structure and becomes smooth. A reasonable estimate 
for the final structure parameter is that of a Plummer 
sphere with a/ ~ 0.75. But the highly fractal initial 
conditions have a structure parameter which is signifi- 
cantly larger than that for a smooth distribution. For 
a fractal with fractal dimension =1.6 numerical experi- 
ments show that ao ~ 1.5 (this is an approximate figure 
and can change dramatically depending on the random 
number seed of the fractal). This means that the degree 
of collapse to virialise from a realistic initial virial ratio 
of Qo = 0.3 is a factor of about 2.5, from Ipc initially to 
around 0.4 pc, which is consistent with our simulations 
(note that the core radius will be significantly smaller 
than 0.4 pc). 

It is important to note that clusters which undergo a 
collapse and dynamically mass segregate will often be sig- 
nificantly dynamically older than their cu rrent crossing 
time might suggest (as also emphasized bv lBastian et al.l 
l2008h . Clusters can dynamically segregate in 1 - 2 Myr 
even if this is comparable to their initial or current cross- 
ing times (as inferred from their sizes) because they have 
undergone a dense phase, which means that they are ac- 
tually many crossing times old (in their cores at least). 

The dense cores are very short-lived as the low-mass 
stars expand due to the transfer of energy from the high- 
mass stars to low-mass stars during the mass segregation. 
The lifetimes of the dense phase is usually only 0.1 - 
0.2 Myr. 

There is one obvious omission from our simulations - 
gas. For simplicity, we neglect the background potential 
of the gas from which the star clusters formed. However, 
we believe that this omission will not affect our results 
significantly. Much of the gas that does not form stars 
has a high velocity dispersion: the cloud as a whole may 
well be unbound due to supersonic turbulence. Stars are 



Adams F. C, Proszkow E. M., Fatuzzo M., Myers P. C, 2006, ApJ, 
641, 504 



able to form cool as they form in converging and shocked 
regions within a globally unbound cloud (e.g. Bate et 
al. 2003; Bonnell et al. 2003). Therefore, we would 
not expect significant amounts of gas to collapse with 
the stars (obviously some will, but the contribution to 
the potential will become small as the stars collapse to a 
de nser c o nfigur ation). This is observed in the simulations 
of iBatd ()2009() . where gas which has not formed stars 
has migrated away from the star forming areas due to 
its initial velocity dispersion. In particular, Fellhauer, 
Wilkinson, & Kroupa (2009) show that the presence of 
a background potential does not significantly alter the 
merging behavior of subclumps. 

4. CONCLUSIONS 

Observations and theory suggest that many young star 
clusters form with a significant amount of substructure. 
Observations also show that young clusters lose their 
substructure on a timescale comparable to their current 
crossing times (1-2 Myr). Simulations suggest that the 
only way in which this could happen is if cl usters are born 
dynamically cool (jGoodwin fc Whitworth 2004). There- 
fore we argue that the correct initial conditions of many 
star clusters must be a cool, clumpy distribution. 

We conduct an ensemble of simulations of cool, clumpy 
clusters (with fractal dimension 1.6 and virial ratio 0.3) 
to investigate the early dynamical evolution of such clus- 
ters. We also find that cool, fractal clusters tend to 
dynamically mass segregate down to a few solar masses 
during a dense, but short-lived, state at the end of their 
collapse. Such limited mass segregation of only the most 
massive stars is what appears to be observed i n the ONC 
(|Allison et al.l[200l iMoeckel fc Bonnel]|l2009D . 

Such rapid dynamical mass segregation in clusters with 
realistic initial conditions shows that the most massive 
stars do not have to form in the centers of clusters for 
very young clusters to be mass segregated. This shows 
that massive stars could form in relative isolation in large 
cores and mass segregate later, possibly avoiding the 
need for competitive accretion to form the most mas - 
sive stars in the center of a cluster (jBonnell et al.l[200l[ l. 
However, competitive accretion may still play an impor- 
tant role, further increasing the masses of the most mas- 
sive stars if they mass segregate while gas is still present. 
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